Lithium Ion Transport Mechanism in Ternary Polymer Electrolyte-Ionic Liquid 
Mixtures — A Molecular Dynamics Simulation Study 
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The lithium transport mechanism in ternary polymer electrolytes, consisting of PE02oLiTFSI 
and various fractions of the ionic liquid PYR13TFSI, is investigated by means of MD simulations. 
This is motivated by recent experimental findings^, which demonstrated that these materials dis- 
play an enhanced lithium mobility relative to their binary counterpart PE02oLiTFSI. In order to 
grasp the underlying microscopic scenario giving rise to these observations, we employ an analytical. 
Rouse-based cation transport model^, which has originally been devised for conventional polymer 
electrolytes. This model describes the cation transport via three different mechanisms, each char- 
acterized by an individual time scale. It turns out that also in the ternary electrolytes essentially 
all lithium ions are coordinated by PEO chains, thus ruling out a transport mechanism facilitated 
by the ionic-liquid molecules only. Rather, the reason for the enhanced lithium mobility can be 
found in the plasticizing effect of the ionic liquid, which enhances the dynamics of the PEO chains 
and thus also the motion of the attached ions. Additional focus is laid on the prediction of lithium 
diffusion coefficients from the simulation data for various chain lengths and the comparison with 
experimental data, thus demonstrating the broad applicability of our approach. 

PACS numbers: 



I. MOTIVATION 

Solid polymer electrolytes (SPEs) are promising candi- 
dates for lithium ion batteries, as they are ideal to create 
small and light-weighted but powerful energy storages^. 
The classical SPEs consist of an amorphous polymer ma- 
trix, e.g. poly(ethylene oxide) (PEO), and a lithium 
salt dissolved in i1>^i^. By using lithium salts with large 
anions such as lithium-bis(trifluoromethane)sulfonimide 
(LiTFSI), the crystallization or salting-out, which is a 
common problem for inorganic salts, can be suppressed 
as the negative charge is delocalized over the whole an- 
ion. However, at ambient temperatures, the conductivity 
of most SPEs is still too low for an efficient technolog- 
ical use. Among several other remedies^"— that have 
been proposed to overcome this deficiency, the incorpo- 
ration of a room temperature ionic liquid (IL) seems to 
be a very promising approachi*ii, as the resulting ternary 
electrolytes show both an increased conductivity and in- 
herent stability. Further benefits directly stem from the 
advantages of the ILs themselves: They are non-volatile, 
non-fiammable^ and exhibit a wide electrochemical sta- 
bility windowi^. 

However, it is not yet fully understood in how far the 
lithium transport mechanism in these materials changes 
relative to the conventional polymer electrolytes. For 
instance, it was speculated^i that the lithium ions be- 
come progressively coordinated by the anions from the 
IL and are thus decoupled from the rather slow PEO 
chains. Alternatively, one might also expect that the IL 
enhances the PEO dynamics and serves as a plasticizer 
in this way, which is a common observation when adding 
low-molecular solvents to PEO-salt systems^"—. In this 



work, we utilize molecular dynamics (MD) simulations to 
unravel which scenario is responsible for the experimen- 
tal observations^. In order to quantify the lithium mo- 
tion, we employ an analytical cation transport mode)^^, 
which has originally been devised for binary polymer elec- 
trolytes. 

In the following, we will briefly review this model. Our 
description is based on both the Rouse model^^ as well as 
the Dynamic Bond Percolation (DBP) modeli^, and dis- 
tinguishes three different microscopic lithium ion trans- 
port mechanisms (see also sketch in Figure [1]): 1. The 
ions diffuse along the PEO backbone to which they are 
attached. This motion can be characterized by the time 
scale Ti the ions need to explore the entire PEO chain. 
2. For ambient temperatures, i.e. T ^ Tg, the PEO 
chains are naturally also subject to thermal motion, car- 
rying the attached ions in this way. In case of Rousean 
motion, the polymer dynamics and thus motion of the at- 
tached ions can be quantified by an effective Rouse time 
T2. 3. Finally, an ion bound to a specific PEO chain 
can be transferred to another chain. The mean residence 
time at a given chain is denoted as T3 in the following. 
As demonstrated by earlier MD studies^, the last mech- 
anism can also be viewed as a renewal process within the 
framework of the DBP modeliS, since the dynamics of a 
given lithium ion becomes independent of its past after 
being transferred to another PEO chain. In the limit of 
long chains, this mechanism therefore also allows the ions 
to cover macroscopic distances. 

Of course, switching to the ternary electrolytes, it 
is a priori unclear if this scenario changes only quan- 
titatively - reflected by different values for ri, T2 and 
T3 - or if the lithium ion transport mechanism also 
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Figure 1: Scheme illustrating the three different cation trans- 
port mechanisms in PEO-salt electrolytes. 

changes on a qualitative level. In particular, we focus on 
two ternary polymer electrolytes with a slightly differ- 
ent IL cation than in the experimentsi, i. e. A'^methyl-iV- 
propylpyrrolidinium (PYR13TFSI), with a stoichiometry 
of PEOzoLiTFSI • 0.66 PYR13TFSI and PEOzoLiTFSI • 
3.24 PYR13TFSI, respectively. The binary polymer elec- 
trolyte, PE02oLiTFSI, serves as a reference substance. 
For convenience, PEO will be abbreviated as 'P' (i. e. 
polymer or PEO) and LiTFSI as 'S' (i.e. salt) in the 
following, leading to the short-hand notation P20S • x IL 
with X — X = 0.66 and x = 3.24. 



II. SIIVIULATION DETAILS 

The simulations were performed with the AMBER 10 
Program Package^. Here, the sander module was 
modified, allowing us to use a many-body polarizable 
force field specifically designed for PEO/LiTFSI^^'^^ and 
PYRisTFSI^". The simulation box contained 10 PEO 
chains with iV = 54 monomers each as well as 27 LiTFSI 
ion pairs, yielding a concentration of ether oxygens (EOs) 
to lithium ions of EO : Li = 20 : 1. Additionally, 
the two ternary systems contained 18 or 87 PYR13TFSI 
molecules, corresponding to x = 0.66 and x — 3.24. After 
equilibration runs of 70—80 ns in the NpT ensemble, pro- 
duction runs with a length of 200 ns have been performed 
in the NVT ensemble, collecting data every picosecond. 
An elementary integration step of 1 fs was used, while the 
systems were coupled to a Berendsen thermostat with 
a reference temperature of 423 K. All bonds involving 
hydrogen were constrained by the SHAKE algorithm^^. 
The induceable point dipoles were integrated by a Car- 
Parrinello-like schemeS^. 



III. STRUCTURAL PROPERTIES 

We find for all electrolytes that virtually all lithium 
ions are coordinated to one or two PEO chains, thereby 
giving a first hint that also for the ternary systems the 
cation transport entirely takes place at the PEO chains. 
The percentages of lithium ions coordinating to two PEO 
chains are 53 % for P20S, 47 % for P2oS-0.66 IL and 24 % 
for P20S • 3.24 IL, thus indicating a dilution effect which 
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Figure 2: MSD of the lithium ions (main panel) and the center 
of mass of the PEO chains (inset). The solid lines in the main 
panel correspond to the model predictions. 

reduces the probability for an ion to encounter a second 
PEO molecule. This is also supported by the observation 
that the amount of Li"*" coordinating two PEO chains 
decreases linearly with x. 

The predominant lithium coordination consists of 4 — 5 
EOs, which is in good agreement with experimental 
datap^i and quantum chemistry calculation o^^'^^ (see Fig- 
ure [5] in the supporting information). In those complexes 
where the 4—5 EOs originate from a single PEO molecule, 
the polymer chain wraps helically around the cation. For 
complexes involving two PEO chains, typically 2 — 3 EOs 
from each chain coordinate to the ion. Additional coordi- 
nations by TFSI oxygens are rather rare (about 12—20 %, 
see Figure , and in most cases the anion coordinates 
only briefiy to the lithium ion. Again, the amount of 
Li^ coordinating to TFSI increases linearly with a;, thus 
indicating that this effect is purely statistical. 

For the conformational properties of the PEO 
molecules, we naturally observe a contraction of the 
chains due to the crown-ether-like coordination sphere 
of the lithium iona^^. This manifests itself by a de- 
crease of the mean squared end-to-end vector (R^) for 
the electrolytes ((R^) = 1662 A^, 1570 and 1571 A^ 
for a; = 0, 0.66 and 3.24, respectively) as compared to the 
pure PEO melt ((R^) = 1979 A^). In case of the ternary 
electrolytes, one might expect that the PEO chains are 
swollen on a global scale, since the addition of IL would 
induce a crossover from a polymer melt to a semidilute 
solution. However, from both (R^) and the scaling of the 
Rouse-mode amplitude a^^'^^ , we observe no such feature. 



IV. DYNAMICS OF THE LITHIUM IONS 

Although from a structural point of view no significant 
differences between the individual electrolytes emerge, 
we observe a clear increase of the overall lithium mean 
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Figure 3: MSD of the average EO index {An^{t)). 

square displacement (MSD), especially for the subdiffu- 
sive regime at t = 1 — 10 ns (crosses in Figure[2]), whereas 
the onset to diffusion essentially occurs on a similar time 
scale, i. e. t = 20 — 50 ns. A similar increase can be 
found for the MSD of the entire PEO chains (inset in 
Figure [2]) . Although the larger center-of-mass motion of 
the PEO chains contributes to the effective transport of 
the lithium ions due to their cooperative diffusion, the 
differences in the inset of Figure [5] can only partly ex- 
plain the increase of the lithium MSD. Thus, at least one 
of the three transport mechanisms must become faster 
due to the presence of the IL, which we will investigate 
in the following. 

A. Motion Along the PEO Backbone 

First, we start with the evaluation of the diffusion 
along the PEO backbone. In order to quantify this 
mechanism, it is convenient to successively number all 
monomers at a given PEO chain in order to express the 
lithium position by the average EO index n, and to cal- 
culate an effective MSD (An^(t)) along this coordinate 
(Figure©. 

Starting from approximately 100 ps, the dynamics is 
slightly subdiffusive, characterized by an exponent of 
a w 0.8. Qualitatively, this behavior is found for all 
IL concentrations as well as for both coordination types 
of the lithium ions (i. e. bound to one or to two PEO 
chains). For sufficiently long chains as well as a signifi- 
cant amount of PEO— Li"'' complexes that exist through- 
out the entire observation time, {An'^{t)) will show dif- 
fusive behavior on longer time scales. Naturally, in the 
limit t — !> oo, one would expect a crossover to a plateau 
for finite chain lengths. From Figure O neither of these 
two effects can be found, indicating that most life times 
of the PEO— Li'*' complexes are too short. Rather, the 
lithium ions only move on average 7 — 8 monomers during 
10 ns. Keeping in mind that the lithium ions are typi- 
cally bound to 4 — 5 monomers, these findings imply that 
the ions have barely left their own coordination sphere 
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Figure 4: MSDs of the average EOs (circles), the bound EOs 
(crosses) and the lithium ions bound to these EOs (inset, solid 
lines). All MSDs have been computed in the center-of-mass 
frame of the PEO chain. The dashed and dotted lines show 
the respective Rouse fits. 

during the accessible time scale. Therefore, finite size 
effects of the PEO chains are irrelevant in the present 
case. 

In total, no significant dependence of (An^(i)) on the 
IL concentration can be observed, indicating that the sur- 
rounding molecules (i. e. PEO chains or IL) have virtu- 
ally no influence on this mechanism. It is interesting 
to note that the magnitude of {An^{t)) is essentially the 
same for lithium ions bound to one or to two PEO chains 
(not shown). 

Although the motion along the chain is still slightly 
subdiffusive, approximate diffusion coefficients Di{t*) 
were extracted from Figure |3] in order to determine an 
estimate for ti. The Di were obtained from {An? it)) 
from the Einstein relation, Di{t*) = (Ar?{t*)} / {2t*), 
where t* represents a specific time scale. Ideally, one 
would choose t* = (determined below) in order to es- 
timate the net effect of this mechanism. Unfortunately, 
the statistics of (An^(t)) are too bad on these time scales. 
Therefore, £)i(t3) was determined from Figure [3] by ex- 
trapolation, and the ri-values in turn were calculated ac- 
cording to^n = (A^-l)V(7r^ Di) (TableU). One observes 
that Ti slightly decreases with increasing IL concentra- 
tion. However, these values also implicitly depend on T3, 
which to a large extend explains the weak dependence on 
the IL content (cf. Figure 121). 



B. Polymer Motion 

Next, we focus on the polymer dynamics itself. Fig- 
ure |4] shows the MSD of the EOs relative to the center 
of mass of the PEO chain. This quantity has been com- 
puted for all EOs (i. e. irrespective of the presence of an 
ion) , for EOs bound to a lithium ion as well as for the re- 
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P20S 


147 


45 


167 


17 


2.5 


3.133 


2.945 


1.947 


P20S • 0.66 IL 


140 


37 


89 


18 


1.0 


3.698 


3.542 


2.309 


P20S • 3.24 IL 


127 


24 


68 


24 


8.5 


4.744 


4.257 


2.392 



Table I: Parameters characterizing the three transport mechanisms (see text for further explanation). 



spective attached ions. The criterion to consider a cation 
or EO as bound was that the average EO index of the 
ion did not change more than one, i. e. | An(i)| < 1 for all 
time frames during t. For the bound EOs, no further dis- 
tinction between additional coordinations of the lithium 
ion to another PEO chain or a TFSI molecule was made. 
Thus, these effects are already implicitly contained in the 
curves in Figure |4l 

The average EOs (circles) show typical Rouse-like mo- 
tion with the characteristic relaxation time tr. The dy- 
namics of the bound EOs (crosses) is qualitatively the 
same but protracted. Therefore, it is possible to char- 
acterize the intramolecular dynamics of the bound EOs 
by a larger, effective Rouse time T2. The lithium ions 
attached to these EOs (inset, solid lines) closely follow 
the bound EOs, which gives clear evidence for their co- 
operative motion. On short time scales, the MSD of the 
EOs is larger than the lithium MSD due to the addi- 
tional internal degrees of freedom of the PEO backbone, 
but the MSD of the bound cations catches up on longer 
time scales (approximately 1 ns). Thus, T2 characterizes 
both the dynamics of the bound PEO segments as well 
as the attached lithium ions. For the ternary systems, 
the respective curves are qualitatively the same as in the 
inset of Figure H] 

Figure H] also shows the Rouse fits, i.e. 
ffR(t) = 2{^l) 7r-2 ^ [1 - exp {-tpVm)] p-' 
for the average (dashed lines) and 32 (i) = 

2(R2) 77-2 ^00 ^ _ gj^p (_^p2/^2)j p-2 foj. i-j^g ^^^^^ 

EOs (dotted lines). Of course, the precise value of tr 
and T2 obtained from these fits also depends on the 
value of (Re) in the prefactor. In order to obtain a fit 
consistent with the plateau value at large t, the MSDs of 
the average EOs were fitted using two parameters, i. e. 
tr and (Rp). Subsequently, the MSDs of the bound EOs 
were fitted using this value in combination with a single 
fit parameter T2 only. The resulting values for (R^) are, 
as already discussed above, approximately constant for 
all electrolytes. The respective tr- and T2-values are 
summarized in Table H] (deviations from our previous 
studyii^ on P20S arise from the shorter simulation length 
of about 27 ns as well as the modified fitting procedure). 

Both Tr and T2 decrease significantly, clearly indicating 
that the dynamics of the PEO segments becomes faster 
with increasing IL concentration. Therefore, the IL can 
be regarded as plasticizer. For the average segments, 
the dynamics for P20S • 3.24 IL is nearly the same as for 
pure PEO (tr = 22 ns), showing that the plasticizing 
approximately cancels with the slowing-down caused by 



the coordinating lithium ions as found for P20S. The 
presence of the IL also enhances the motion of the bound 
segments, and, as a result, the dynamics of the respective 
attached lithium ions, leading to an increase of the overall 
lithium MSD. Here, experimental studies reported simi- 
lar findings for other plasticizers like ethylene/propylene 
carbonate^iiS or short PEO chains embedded in a high- 
molecular weight matrix^jlS. 

Finally, the center-of-mass dynamics of the PEO chains 
becomes also faster with increasing IL concentration (in- 
set of Figure [2]). Thus, for finite chain lengths, the plasti- 
cizing effect is twofold since both the internal, segmental 
PEO dynamics and the center-of-mass motion is accel- 
erated, which both contributes to the lithium diffusivity 
(see discussion below). 

It should be pointed out that it is questionable if 
cations bound to two PEO chains show the same dy- 
namical behavior as ions bound to one chain only, since 
the former complexes could be regarded as transient 
crosslinks, which would significantly impede the polymer 
dynamics. A more detailed analysis (see supporting in- 
formation) indeed revealed that there is a fundamental 
difference between these two coordination types. While 
for lithium ions coordinated to a single PEO chain ions 
that diffuse along the backbone are significantly faster 
than those that remain bound to the same monomers, no 
difference between these two categories could be found for 
cations attached to two PEO chains. This implies that 
the cations bound to two PEO chains experience no ef- 
fective transport due to the diffusion along the chain. 
Rather, the PEO chain moves reptation-like along its 
own contour past the cation, which results in a non- 
zero (An2(t)) (Figure [3]), but does not contribute to the 
lithium transport. However, this feature can be easily 
implemented by weighting the effect of ri according to 
the ratio of ions bound to one PEO chain only (see dis- 
cussion in the supporting information). 



C. Interchain Transfer 

Next, we will determine the renewal time for the in- 
dividual electrolytes. The number of transfer processes 
A'tr was counted from the simulations, and the Ta-values 
were determined according to T3 = imaxA^Li+ /A'tr, where 
imax = 200 ns is the simulation length and Ni^i+ — 27 
is the number of lithium ions in the simulation box. Of 
course, it is questionable if brief transfers followed by suc- 
cessive backjumps to the previous polymer chain serve 
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as a renewal process in the strict sense, since the hthium 
dynamics will not become uncorrelated to its past after 
such an event. A more detailed analysis (not shown) re- 
vealed that these non-Markovian, short-time backjumps 
occurred up to 100 ps, which we used subsequently as a 
criterion to define real renewal events. In cases where the 
transfer was mediated by TFSI anions only (probability 
piL in Table H]) , we found that the displacement the ion 
covers in the IL-rich region was sufficiently small, so that 
the contribution of these transfers to the lithium MSD is 
negligible. 

The individual values for T3 are summarized in Table ID 
One observes that T3 increases with increasing IL concen- 
tration. Since the PEO molecules become more and more 
diluted, this can mainly be explained as a concentration 
effect. Obviously, the critical step for a transfer process 
is the encounter of a another PEO segment. 



V. APPLICATION OF THE TRANSPORT 
MODEL 

As a first consistency check of our description, we 
employ the underlying physical picture of the transport 
model to reproduce the lithium MSD in Figure [21 During 
the residence time i at a given PEO chain, the MSD 512 
of the lithium ion is given by a Rouse-like expression^ 
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where t{^^ = Tj"^ -|- r^"'^ is a combined relaxation rate 
due to both intramolecular mechanisms. After a renewal 
process (i.e. an interchain transfer), the ion dynamics 
becomes uncorrelated to its paslj^ii^, and the motion at 
the new chain is again characterized by Eq. [TJ Thus, 
it is possible to interpret the overall lithium dynamics 
as a random walk, in which the elementary step length is 
given by Eq.[T] and the direction of motion for the succes- 
sive step changes randomly after each renewal event. In 
the following, we assumed that the number of ion trans- 
fer processes during time interval t is given by a Poisson 
distribution with mean t/r^, leading to exponentially dis- 
tributed i at the individual PEO chains. For a given ob- 
servation time t, the lithium MSD due to the three trans- 
port mechanisms was then obtained from the numerical 
average over a large number of Poisson processes. The 
third ingredient required for the total lithium MSD as 
depicted in Figure [2] is the center-of-mass motion of the 
PEO chains, which was directly extracted from the simu- 
lations (inset of Figure [2]) and added to the model curve. 
The resulting predictions are shown as solid lines in the 
main panel of Figure [2] 

For P20S, one observes a nearly perfect agreement 
throughout the entire observation time, starting from 
the Rouse regime to the diffusive long-time limit. This 
demonstrates that the simplistic picture of our model in- 
deed captures the underlying, much more complex micro- 



scopic scenario. Interestingly, essentially the same model 
curve is obtained when uniformly distributed t are used. 
Apparently, a renewal event is unlikely on short time 
scales, and the MSD reduces to Eq. [TJ whereas for large 
time scales only the average residence time T3 is impor- 
tant. 

In case of P20S • 0.66 IL and P20S • 3.24 IL, the model 
curves agree with the empirical lithium MSD for time 
scales larger than about 1 — 2 ns. Slight deviations can 
be attributed to the large uncertainties of the MSD of 
the PEO chains. However, the model prediction system- 
atically overestimates the MSDs of P20S • 0.66 IL and 
P20S • 3.24 IL for short time scales. Here, a more detailed 
analysis (to be published under separate cover) revealed 
that these deviations are caused by hydrodynamic inter- 
actions arising from the presence of the IL. However, on 
larger length and time scales, these hydrodynamic inter- 
actions are screened, which has also been reported for 
other semidilute polymer solutions^^-. Thus, both the 
Rouse-like behavior on intermediate time scales and the 
diffusive regime are correctly reproduced, which clearly 
demonstrates the applicability of our model to the exper- 
imentally relevant long-time limit. 

Motivated by these findings, we finally use the trans- 
port model to compute the total lithium diffusion coef- 
ficient Dli- One of the main benefits of our model is 
the fact that due to its analytical nature Du can be cal- 
culated for all chain lengths N, which immediately fol- 
lows from the scaling laws of the three time scales^, i. e. 
Ti cx N^, T2 cx N"^ and T3 oc N°. 

Of course, for the scaling of T2, entanglement effects 
may become relevant, which would slow down the seg- 
mental dynamics. However, if T3 < Te, meaning that the 
lithium ion leaves the PEO chain before the latter begins 
to reptate (the onset of this regime being characterized 
by the entanglement time Te), the overall dynamics is still 
Rousean^®. For such a scenario, our model can still be 
used to calculate Dli- In the case of PEO, experiments 
revealed that the entanglement regime sets in from about 
N 75^1^. Based on these observations, one can es- 
timate To according to Tc = tu,{N = 75) = rR(75/54)2. 
For P20S, this leads to w 87 ns, which is substantially 
larger than T3. Also in case of the highly plasticized 
P20S • 3.24 IL one finds « 46 ns > t^. Therefore, 
the lithium ion leaves the PEO chain before the tube 
constraints become noticeable, and the formalism of the 
transport model can thus also be applied in the limit of 
long chains. 

As an alternative to the numerical evaluation of the 
Poisson process (model curves in Figure [2]), it is also pos- 
sible to derive an analytical expression for the diffusive 
long-time limit, characterized by Dli. From the distri- 
bution of the residence times i at the individual PEO 
chains, i. e. p{i) = t^^ exp(— i/r^), the mean elementary 
step length {gi2{T3, N))m3, averaged over all possible i, 
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can be obtained: 

(512(7-3, A^))m3 = — / di exp[ — -] gi2{i) 
Jo \ T3/ 
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Since ri and r2 depend on the chain length, the av- 
eraged MSD between two successive renewal events, 
(512(^3, A^))m3, is also a function of N. Moreover, one 
has to take into account the additive contribution due 
to the center-of-mass diffusion of the PEO chains. By 
making use of the renewal property, the overall lithium 
diffusion coefficient can thus be expressed as 



Du{N) 



(512(7-3, ^))m3 

6t3 



DpEoiN), 



(3) 



where I?peo(-^) is the diffusion coefficient of the PEO 
chains with length N. By solving Eq. [5] numerically and 
inserting the result into Eq. [Sj values for Du were ob- 
tained for = 54 {D^=^^) and iV ^ 00 {D(^--^°°), which 
are summarized in Table HI 

For we observe that the values obtained from 

Eq. |3] are in reasonable agreement with the respective 
diffusion coefficients Z?^^^ directly determined from Fig- 
ure [5J again showing that the predictions of our model 
are consistent with the numerical data. When comparing 
D^r^°° with the PFG-NMR datai (D^^ = 0.052 A^ns"! 
for P20S • 1.0 IL, 0.118 A^ns-i for P20S • 2.0 IL and 
0.126 A^ns-i for P20S • 4.0 IL at T = 323 K), one ob- 
serves identical trends as in Table HI namely a significant 
increase of Dli with increasing IL amount. As discussed 
above, these findings can be attributed to the plasticizing 
effect of the IL. 

It is illuminating to look at these values in more de- 
tail, since the segmental motion entering (512(7-3, iV))M3, 
the renewal time T3 and the center-of-mass motion, mea- 
sured by -DpEO, alters the value of Du (Eq.[2]). For these 
reasons, it is desirable to quantify in how far these three 
trends contribute to the overall change of D^. In case of 
iV — ;> 00, the intramolecular contribution (5i2(t3, A^))m3 
increases by about 28 % for P20S • 0.66 IL and 73 % for 
P20S • 3.24 IL, mainly as a result of the increased seg- 
mental mobility (Figure 2]). In contrast to this, the con- 
tribution of the renewal processes decreases by about 
7 % and 29 %, respectively. Nevertheless, the enhance- 
ment of (512(7-3, iV))M3 overcompensates this effect, and 
the overall lithium diffusivity increases by about 19 % 
for X = 0.66 and 23 % for x = 3.24. On the con- 
trary, the changes in T3 dominate in the short-chain 
limit, since the differences in ti and T2 for the individ- 
ual electrolytes disappear for iV ^ 0. Of course, this 
does not imply that also the differences in Du disap- 
pear, as the center-of-mass diffusion of the PEO chains 
becomes faster with decreasing N and with increasing IL 
concentration (Figure [5]). Thus, the stronger increase 
of D^=^^ with X (20 % for x = 0.66 and 45 % for 



X = 3.24, Table H]) as opposed to D^~^°° arises to a sub- 
stantial degree from the contribution of DpEO- In fact, 
for N = 54, (512(7-3, iV))M3/(6r3) increases only by 13 % 
for P20S • 0.66 IL and by 9 % for P20S • 3.24 IL, whereas 
-DpEo becomes faster by about 30 % and 92 %, respec- 
tively. 

So far, we focused on the high-temperature limit which 
we can address in our simulations. Interestingly, it was 
observed that in the low-temperature regime the rela- 
tive increase of D'^^ upon the addition of IL becomes 
even more pronounced^ (i. e. D'j^^ increases by a factor 
of 2.4 when going from P20S to P20S ■ 4.0 IL at 323 Kl, 
whereas the relative increase of D^f~^°° between P20S and 
P20S • 3.24 IL is only haff as large (Table HI). Affhough 
simulations at low temperatures would be too costly, one 
would nevertheless expect that the plasticizing effect at 
least partly accounts for the larger relative increase of 
Du in this regime. This is also indicated by DSC mea- 
surements^, which revealed that the glass-transition tem- 
perature Tg decreases significantly (up to 35 K) with the 
IL fraction. Thus, it is likely that the enhanced poly- 
mer dynamics, apart from other, possibly related effects, 
plays a crucial role at low temperatures. 



VI. CONCLUSION 



In this study, we have examined the microscopic 
lithium ion transport mechanism in ternary polymer elec- 
trolytes consisting of PE02oLiTFSI and PYR13TFSI. In 
particular, we addressed the question in how far the mi- 
croscopic scenario of the ion transport changes upon the 
addition of IL, and how the experimentally observed in- 
crease in the lithium ion diffusion coefficient^ can be 
understood in terms of the individual transport mech- 
anisms. To this purpose, an analytical cation trans- 
port model^ was successfully applied. It turned out that 
virtually all lithium ions were coordinated to the PEO 
chains, thus ruling out a transport mechanism in which 
the lithium transportation is decoupled from the poly- 
mer chains. Rather, the main reason for the increase of 
the lithium diffusion coefficient is the plasticizing effect of 
the IL, which enhances the segmental motion of the PEO 
chains, and, consequently, the dynamics of the respective 
attached ions. Another minor counteracting effect was 
the successive dilution of the electrolyte due to the IL, 
which slightly decreases the rate of interchain transfers. 
In the sum however, the plasticizing is dominant, and 
the overall lithium diffusivity increases. For the design 
of novel battery materials, our findings therefore imply 
that a polymer electrolyte which is both highly plasti- 
cized and exhibits a high transfer rate, e. g. facilitated 
by a more coordinating IL, would yield optimal results. 
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^ — PEO^gLiTFSI 

I ' ^ — PEOjgLiTFSI 0.66 IL 

\ _ PEO"xiTFSI 3.24 IL 




Figure 5: Probability p{n) to find a certain coordination num- 
ber n of EOs (irrespective if tlie ion is tied to one or two PEO 
cliains) or TFSI oxygens. 



0-OLiat 1 PEO (M2 only) 




f[ps] 

Figure 6: MSDs of the lithium ions bound to one or two 
PEO chains in PE02oLiTFSI. A second distinction was made 
between cations that diffused along the PEO chain (Ml and 
M2) and ions that remained bound to the same EOs (M2 
only). All curves have been computed in the center-of-mass 
frame. For the other electrolytes, the scenario is qualitatively 
the same. 



VII. SUPPORTING INFORMATION 
A. Coordination Environment 

Figure [5] shows the probability distribution functions 
p{n) to find a lithium ion with n EOs or TFSf oxygens 
in its first coordination shell. 



B. Influence of the Lithium Ion Coordination on 
the Segmental Motion 

In order to elucidate in how far the lithium ions coor- 
dinated to two PEO chains act as temporary crosslinks, 



system 1 PEO [%] 2 PEO [%] 



P20S 




47.2 


52.7 


P20S 


0.66 IL 


53.0 


47.0 


P20S 


3.24 IL 


75.8 


24.0 



Table II: Ratio of lithium ions coordinating to one or two 
PEO chains. 

Figure |6] shows the lithium MSD of lithium ions bound to 
one and to two PEO chains for P20S. A second distinc- 
tion was made if the ions diffused along the PEO chain 
(both intramolecular mechanisms) or remained bound to 
the same EOs (polymer dynamics only, here, the same 
criterion as for Figure IH \An{t)\ < 1, has been ap- 
plied). In case of cations bound to one PEO chain only 
(Figure E]), ions undergoing both types of intramolecular 
transport are clearly faster than those remaining close 
to the initial position on the PEO chain. For lithium 
ions connected to two PEO chains, no significant dif- 
ference between these two categories can be observed 
in the MSD (Figure |6]). This implies that the cations 
bound to two PEO chains experience no effective trans- 
port due to the diffusion along the chain. Rather, the 
PEO chain moves reptation-like along its own contour 
past the cation, which results in a non-zero (An^(i)), 
but does not contribute to the overall lithium transport. 
For all other investigated electrolytes, the observations 
from Figure ini are qualitatively the same. 

With respect to the transport model, in particular 
Eq. [1] this effect can be easily captured. Here, the only 
additionally required parameter is the ratio of lithium 
ions coordinated to one or two PEO chains (Table lll|): 

gi2 {t) = ri PEO 512 (t) + - ri PEo) 92 (t) (4) 

For the fraction riPEo of cations bound to one PEO 
chain, Eq. [1] remains valid, whereas for ions bound to 
two chains with ti^2PEO — ^ 00 only T2 is important. 
Within this approximation, the average dynamics of the 
lithium ions (i. e. averaged over ions bound to one or two 
PEO chains while simultaneously allowing the interme- 
diate exchange between these two coordination states) is 
estimated from the structural property ri peo only. 

Figure [7] shows that this rather simplistic picture is 
indeed valid to a good approximation. Here, Figure [7] 
displays the MSD of ions bound to one or two PEO 
molecules in P20S as extracted from the simulations. 
These curves have in turn been used to calculate an ap- 
proximate average MSD according to Eq. 01 which is also 
shown in Figure [71 In fact, the agreement of gi2{t) with 
the average MSD (i. e. a lithium ion that remained on the 
same chain irrespective of other coordinations) is nearly 
quantitative. Deviations for larger time scales are due 
to bad statistics. In the same spirit, the inset of Fig- 
ure [7] shows the ratio between the approximate gi2{t) 
and the average MSD directly calculated from the simu- 
lations for all systems. As for P20S, this ratio is close to 
unity for all other electrolytes. These observations high- 
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Figure 7: MSD of lithium ions bound to one or two PEO 
chains as weU as the MSD of cations which are bound to the 
same PEO molecule during t (irrespective of further coordi- 
nations) for P20S. The MSDs of ions bound to one and to 
two chains have been used to compute an approximate aver- 
age gi2{t) according to Eq. U) Inset; ratio of 512 (i) relative to 
the real average MSD as extracted from the simulations. 



light that, since the crosslinks are temporary and an ex- 
change between both coordination types takes place, the 
long-time behavior of the intramolecular dynamics may 
be estimated by the average in Eq. 31 as nicely demon- 
strated in Figure [71 



